Modeling the Antiferromagnetic Phase in Iron Pnictides: Weakly Ordered State 
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We examine electronic states of antiferromagnetic phase in iron pnictides by mean-field calcu- 
lations of the optical conductivity. We find that a five-band model exhibiting a small magnetic 
moment, inconsistent with the first-principles calculations, reproduces well the excitation spectra 
characterized by a multi-peak structure emerging below the Neel temperature at low energy, together 
with an almost temperature-independent structure at high energy. Investigating the interlayer mag- 
netoresistance for this model, we also predict its characteristic field dependence reflecting the Fermi 
surface. 

PACS numbers: 75.30.Fv, 78.20.Bh, 72.80.Ga 
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Introduction. — Iron pnictide superconductors have re- 
ceived considerable attention since the discovery of super- 
conductivity in LaFeAsOi-^F^ [l[. The superconductiv- 
ity appears next to the antiferromagnetic (AFM) phase 
that emerges when x = for LaFeAsOi-^Fa; with the 
Neel temperature Tjy ~ 137 K Since the pairing of 
electrons is expected to be mediated by magnetic fluc- 
tuations [1, H[ , the understanding of the AFM phase is 
crucial for clarifying the mechanism of the superconduc- 
tivity in the iron pnictides. 

First-principles calculations of the band structure are 
considered to have enriched our knowledge of the un- 
derlying electronic structures in the iron pnictides [5j. 
Those calculations, however, disagree with experiments 
on the magnitude of the magnetic moment in the AFM 
phase @, Q: The calculations yield a strongly ordered 
AFM state with the magnetic moment of ~ 2 /iB, which 
is larger than the experimental ones of < 1 [Xb [2j. To 
this discrepancy in magnetic moment, several interpreta- 
tions have been provided, such as the presence of domain 
motions dynamical spin fluctuations and the re- 
quirement of negative U [9j. Still unknown is whether 
the weak order model is appropriate — though we know 
the strong order model seems inappropriate for the iron 
pnictides — since the properties of weakly ordered states, 
apart from the small magnetic moment, have not been 
much studied yet. 

The magnetic order strength affects not only the mag- 
netic excitations but also the charge excitations because 
of the presence of gap (not a full but a partial gap). 
In fact, the in-plane optical conductivity for BaFe2As2 
and SrFe2As2 as well as EuFe2As2 [ll[ has shown 
significant change of its spectra with reducing temper- 
ature across Tjy. A new multi-peak structure appears 
at around 0.05-0.18 eV depending on material, while the 
spectra above ~ 0.6 eV is almost unchanged. A similar 
peak structure has also been observed in polycrystallme 
LaFeAsO [HJ . It has not yet been clarified whether such 
a spectral change is consistent with the emergence of the 
weak magnetic order. This is crucial for modeling the 
AFM phase as well as the superconducting phases of the 
iron pnictides. 



In this Letter, we discuss the excitation properties of 
the weak order model. Investigating the optical conduc- 
tivity arising from interband transitions for weakly or- 
dered states, we verify that the weak order model re- 
produces the experiments and depict the system well. In 
addition, we predict the field dependence of the interlayer 
magnetoresistance in this system. 

Mean-field five-band model. — Considering an Fe 
square lattice, we start with the Hamiltonian for a d- 
electron system H = Hq + Hj. Here, 



Hn = 



E 



t(A x ,A y ;^,v)e 



ik-A 



is the five-band hopping Hamiltonian, where c^ vcr cre- 
ates an electron with a wave vector k and a spin a at 
an orbital a, the hopping energy t(A x , A y ; fj,, v) is given 
from Ref. Q, and A = (A x ,A y ). Hi is the interaction 
Hamiltonian (l3j : 

nV > 2U-5J 

ni = U y ^Tli^nial H 2^ ni^Uiuai 



(2) 



where i is the Fe-site index; U, the intra-orbital Coulomb 
interaction; J, the Hund coupling. The number and the 
spin operators are defined as = ct and Si M = 

\ Stt' c lfj,T^TT'Cifj,T' , respectively, with a the Pauli spin 
matrix. We set the Fe-Fe bond length a to be unity and 
x and y to be along the nearest Fe-Fe bond directions. 

To obtain a model of the ordered state in iron pnic- 
tides, we solve mean-field equations selfconsistently in 
the same manner as in Ref. [14]. Considering the mag- 
netic ordered states with the spin-density-wave (SDW) 
ordering vector Q = (tt, 0), we take the order parame- 
ters such as (n Qp „ CT ) = i Dk( c k+Q M( r c k^a} with N the 
number of k points in the first Brillouin zone of the five- 
band paramagnetic system. We obtain the quasiparticlc 

state "fLa = E;=o,iE /1 Vw( k + i Q)4+/Q i u CT with the 
energy E ke!T . 
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FIG. 1. (Color online) (a) The magnetic moment per Fe site 
for different parameter sets: U = 1.0 eV (squares), 1.1 eV 
(circles), 1.2 eV (triangles), and 1.4 eV (diamonds). Fermi 
surfaces are plotted for the weak order case (b) and the strong 
order case (c) in the magnetic Brillouin zone of the five-band 
model (-| < < |, -1 < < 1). The labeled Fermi 
surfaces (a, /3, and 7) contribute to the interlayer magnetore- 
sistance (see text). 



In our calculation, we change the parameters U and J 
to control the order strength expressed by the magnetic 
moment M = ^("Q^t ""Qwll/'B with \x B the Bohr 
magneton. The calculated magnetic mometns for several 
parameter sets are shown in Fig. [TJa) . The computations 
are performed on a system with TV = 400 x 400. 

For the subsequent arguments, we take the following 
two parameter sets, which give the magnetic moments 
M = 0.4^b and 2.3/zb [Fig. Q] (a)], as models of weakly 
and strongly ordered states, respectively: U = 1.1 eV, 
J = 0.2 eV; and U = 1.4 eV, J = 0.3 eV. The former 
agrees well with the experiments @ on the magnetic mo- 
ment. The latter agrees well with the first-principles cal- 
culations not only on the magnetic moment but also on 
the Fermi surface and the partial density of states [ll|. 
For the parameter sets that give similar magnetic mo- 
ments, we obtain results with similar features to those 
discussed in this Letter. 

The magnetic ordering is achieved with a SDW gap; 
this is not a full gap, and the system stays metallic. The 
Fermi surfaces are drawn in Figs.[ljb) and (c), the data of 
which are used to calculate the interlayer magnetoresis- 
tance below. For other parameter sets, the states with a 
similar magnetic moment possess a similar Fermi surface. 
In all cases, the circular Fermi surface at (0, 0) is hole- 
like and others are electronlike. There are two circular 
Fermi surfaces around (0,0) in the paramagnetic case [4(; 
the SDW gap partially opens along the outer (electron- 
like) one owing to the Fermi surface nesting, and small 
segments remain at a or 7. 

In the weak order model, the states on the hole pockets 
at (0, 0) have more the d yz and the d zx characters than 
the other orbitals; especially, the d yz (d zx ) character is 
dominant near k x — (k y = 0). In the electron pocket 
a in Fig. [Tfb), the d zx character is dominant along the 
flat region on the side near (0,0); the d yz character is 



dominant, on the other side. Most of the states on f3 in 
Fig. [ljb) have the d xy character, and the other states on 
(3 show the d xy and d zx characters. The Fermi surfaces 
observed in a recent angle-resolved photocmission spec- 
troscopy study on BaFe2As2 [l|| are similar to a and 
(3- 

In the strong order case, the hole pocket at (0,0) has 
the d zx character dominantly all around it. The states 
on 7 in Fig. [Tf c) possess a mixed character of several 
orbitals: The states on the (0,0) side show rather the 
d x 2_ y 2 character; on the other side, the d zx and d xy char- 
acters. 

Optical conductivity. — The optical conductivity ex- 
periments on iron pnictides have revealed the excitation 
structure in the magnetically ordered state, whose char- 
acteristics become clear by comparing its spectra with 
those in the paramagnetic state at T = 300 K. The para- 
magnetic state shows a broad peak at around 0.5 eV, 
and no peaks below this apart from the Drude peak [Io| . 
Compared with the paramagnetic state, the magneti- 
cally ordered state shows the following features [10| : (i) 
an emergence of the excitations related to the SDW 
gap at around 0.1 eV, (ii) a slight energy shift of the 
0.5 eV structure, and (iii) almost the same structure 
above 0.5 eV. We use these features (i)-(iii) as criteria 
of judgement in evaluating the model of the AFM state 
in Fc pnictide system. 

Now we calculate the optical conductivity for the weak 
order model to examine whether the weakly ordered state 
can model the AFM state. We also discuss the strong 
order model from the same view point. The real part of 
the optical conductivity is obtained as follows: 



Re{er Q/3 (w)} = 



ANuj 



er,k,e,e' 



<d[d]'^-^-"). ( 3 ) 



where / is the Fermi distribution function, e is the ele- 
mentary charge, and Ck«v ^ s defined as 



Ac) _ \p 



A (a) t(A;M,^)( E + W£ev(k + *Q) 

(=0,1 



{(1 - s„ v ) cos[(k + ZQ) • A] - i(l + s„ v ) sin[(k + ZQ) • A] } 

(4) 

with A( Q ) the a component of the vector A and (= 
±1) defined as t(Aji;fi, v) — s (tt „t(Ay; fx, v). We focus 
on the interband transitions but no Drude component. 

The calculated optical conductivities for the direction 
of a = f3 = (1, 1) are plotted in Figs. Ufa)- (c). Within our 
calculations, the temperature effect is only to broaden the 
excitation peaks through the Fermi distribution function 
/, so that we set T = and neglect the effect in the 
following calculations. (We do not argue the experimen- 
tally observed gap-like behavior above Tm below around 
0.6 eV, the so-called psuedogap feature [ToMT^ . fl^ : this 
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FIG. 2. (Color online) Optical conductivity and density of 
states for the paramagnetic [(a) and (d)], the weakly ordered 
[(b) and (e)], and the strongly ordered [(c) and (f)] cases. The 
inset in (a) shows the whole excitation spectra of interband 
transitions. Some characteristic excitation structures and a 
gap are labeled such as Al, A2, etc. The optical conductivity 
is normalized by . The insets in (d)-(f) show the whole 

structure. 



is possibly due to the spin-fluctuation effects [18J beyond 
our mean-field calculation.) 

First, we discuss the spectra in the paramagnetic case 
[Fig.[^a)]. There are two excitation structures: One lies 
in the energy range 0.1-0.6 eV (labeled as Al); the other, 
above 1.0 eV (A2 in the inset). The double broad-peak 
structure below and above 1.0 eV is also observed in the 
experiments on arsenic compounds [l(| and phosphorus 
ones [l9j]. However, the Al and A2 have a fine struc- 
ture [Fig. Ufa)]; the appearance of the structure differs 
from the single-broad-peak structure in the experiments 
around 0.5 eV. In addition, the intensity ratio of A2 to 
Al is much higher than that in the experiments. These 
differences are associated with effects not included in our 
calculation, such as spin-fluctuation effects. Except for 
these details, the broad aspect seems consistent with the 
experiments. 

Next, we argue the weak order case [Fig. HJJb)], com- 
paring with the paramagnetic one. In the weakly or- 
dered case, the excitations related to the SDW gap ap- 
pear around 0.1 eV (Bl). The other parts of the exci- 
tation spectra are almost the same as that in the para- 
magnetic case, except for the small shift from Al to B2. 
These features agree well with the experimental results 
on the criteria (i)-(iii). It is natural that these three fea- 
tures appear because the weakly ordered state has an 
electronic structure similar to that of the paramagnetic 
one [see Figs. [5Jd) and (e)], and the difference lies only 
near the Fermi level, where the SDW partial gap opens 
[indicated as El in Figs.[2fe)]. The interband transitions 



that contribute the Bl excitations occur at the SDW gap 
near (tt/4, 0), which mainly involve the orbital transition 
from d xy to d yz . 

In the strongly ordered case, on the other hand, the 
density of states and the excitation spectra differ from 
the paramagnetic case [see Figs, [^c) and (f)]. The ex- 
citations related to the SDW gap appear around 0.2 eV 
[CI in Fig. IDJc)], whose energy is inconsistent with the 
criterion (i). Other excitation peaks lie above 0.5 eV, 
whose energies appear to be shifted up by > 0.5 eV as 
compared with Al and A2 (the energy region > 1.0 eV 
not shown here) . This disagrees with the criteria (ii) and 
(iii). Therefore, as is already suggested from the mag- 
netic moment study, the strong order model is certainly 
not a valid model in terms of the excitation properties as 
well. 

Here, we mention three dimensional effects. The first- 
principles calculations show three dimensional Fermi sur- 
faces especially in the systems such as BaFe2As2 and sug- 
gest those effects be important. We confirmed that the 
above behavior of the optical conductivity, however, is 
unchanged even in the presence of the d-d hopping in the 
z direction (not shown here). 

Hence, our optical conductivity calculations suggest 
that the weakly ordered state be a valid model of the 
ordered state in iron pnictides. 

Interlayer magnetoresistance. — For a further experi- 
ment to examine the weak order model, we propose an 
interlayer magnetoresistance experiment. Below, calcu- 
lating the interlayer magnetoresistance, we provide an 
aspect of its magnetic-field and azimuthal- angle depen- 
dence for the benefit of experimental research. The in- 
terlayer resistance at zero temperature p Z z{4>) under a 
magnetic field B = B(cos(^>), sin(0), 0) parallel to the 
layer is given by [20l ] 

i _e^_ / tcoA N sr^ f rflk r 

with t c the interlayer tunneling, c the speed of light 
in vacuum, and a c the interlayer distance, where 3 = 
(vketr xB) + T 2 , and impurity scattering effects 
upon interlayer hoppings are included by the constant 
parameter T = H/(2t) with r the scattering time. Here- 
after, we set r = 2.0 x 10~ 4 eV. 

In Fig. [3J the calculated interlayer magnetoresistances 
are plotted. In the weak order case [Fig. [5Ja)], as B 
increasing from 2 to 6 T, the local minimal value of 
Pzz(4>) I Pzz(0) decreases first and then increases. Further- 
more, at 8 T, p Z z{4>)l Pzz(0) shows a double-peak struc- 
ture. 

Here, we interpret the formula ([5]) to explain the above 
features of p zz (<j>)/p zz (0). It follows from the formula ([5]) 
that a Fermi surface with a flat region leads to an oscil- 
lating behavior in the cf> dependence of p zz ((f)) / p zz (0) [2(| • 
Contrary, a circular Fermi surface hardly contributes to 
the <f> dependence. In fact, the <f> dependence mainly 
arises from the two Fermi surfaces a and (3 in Fig. [TJb) 
(and their symmetric counterparts); in addition, each 
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FIG. 3. (Color online) The azimuthal-angle dependence of 
the interlayer magnetoresistances for the weak order model 
(a) and the strong order (b) under magnetic fields B — 2 T 
(solid), 4 T (dotted), 6 T (dashed), and 8 T (dash-dotted). 



shows a different <f> dependence as well as a different B 
dependence. 

For a small B (<C 2 T), the formula © is almost inde- 
pendent of <p (note that S _1 w T 2 ), and the <j> dependence 
of p zz is negligible. The <p dependence emerges when B 
becomes larger than a threshold, which depends on the 
Fermi velocities and T — the larger the Fermi velocities, 
the lower the threshold. By investigating the Fermi ve- 
locity, we find the threshold for a is lower than (3. 

For B = 2 T and 4 T, the a contribution is domi- 
nant and tends to generate a valley in p zz {(f)) / p zz {Q) at 
4> = 90, owing to the flat segment of a. The (3 contribu- 
tion, on the other hand, emerges for B > 6 T and tends 
to develop a hump at <j> = 90 owing to the flat segment 
at around (n, 0.27r); this suppresses the valley structure 
there. Hence, the crossover of the dominant contribu- 
tion of the Fermi surface results in the non-monotonic B 
dependence. Also, the double-peak structure at 8 T in 
Fig. [21a) is attributed to the total contribution of a and 
(3. For larger B, the hump component due to (3 at <f> — 90 
develops more since the flat segment of (3 is longer than 
that of a. 

For the strong order model [Fig.[3fb)], in contrast, the 
4> dependence at B — 2 T is 90-degrec shifted from that 
for the weak order model. In addition, we find a clear 
difference in the B dependence: p zz ((j>)/ p zz (0) changes 
monotonically unlike the weak order case. This is be- 
cause only one type of Fermi surface — 7 in Fig. QJc) — 
contributes to it in the strong order model. 

Hence, we predict a non-monotonic B dependence of 
the interlayer magnetoresistance reflecting the presence 
of the Fermi surface (3 in the weak order model, rather 



than a monotonic one in the strong order model. We add 
that, with a sufficiently weak scattering, the double-peak 
structure at ~ 8 T would be observable — this structure 
may appear as a dip structure for larger B. These two 
features (the non-monotonic behavior and the double- 
peak structure) are the crucial consequence of the multi 
Fermi surfaces (a and (3). Although the above argument 
is based on the calculations with the scattering rate that 
is the same for each band and isotropic, we can discuss 
it from the same point of view even for the case of dif- 
ferent scattering rates. From calculations with different 
scattering rates for each band (not shown here), we veri- 
fied that at least one of the two features arising from the 
multi Fermi surface appears in the weak order case. In 
addition, we examined the case of the anisotropic scatter- 
ing rate [2l[. The anisotropy causes ripples on the angle 
dependence of the interlayer magnetoresistance; still, the 
features arising from the multi Fermi surface is robust 
against the presence of the anisotropy. 

Conclusions. — We have investigated the weakly or- 
dered state as a model of the striped AFM state in the 
parent compounds of iron-based superconductors. Com- 
pared with the experiments of the optical conductivity, 
our calculations suggest that the weak order model well 
reproduce the characteristics of the experimental results 
in terms of the criteria (i)-(iii). We have also verified that 
the strong order model, which corresponds to the model 
obtained from the first-principles calculations, does not 
reproduce the experiments. From these calculations, we 
conclude that the weak order model is a valid model of 
this system. In addition, we have calculated the inter- 
layer magnetoresistance for the models. From the cal- 
culation results, we predict that the interlayer resistance 
changes non-monotonically as the magnetic field increas- 
ing up to 8 T and a double-peak (or a dip) structure un- 
der large magnetic fields. Further experiments on inter- 
layer magnetoresistance are desired to validate our con- 
clusion of the weak order model. 
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